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We numerically simulate the dynamics of aggregation of interacting atomic clusters deposited on 
a surface. We show that the shape of the structures resulting from their aggregation-limited random 
walk is affected by the presence of a binary interparticle Laplacian potential due to, for instance, the 
surface stress field. We characterize the morphologies we obtain by their Hausdorff fractal dimension 
as well as the so-called external fractal dimension, which appears more sensitive to the potential. 
We demonstrate the relevance of our model by comparing it to previously published experimental 
results for antymony and silver clusters deposited onto graphite surface. 



I. FREE AND BIASED DIFFUSION MODELS 

Highly ramified dendritic structures often appear in 
nature in such different situations as cell colony growth 
dielectric breakdown [2], viscous fingering in Hele- 
Show cells [3|, liquid crystal growth Q, electro-deposition 
(Ell , single-atom [6|, @] and atomic cluster depositions 



l!fl. Though the objects observed are not mathemati- 
cally rigorous fractals, the self-similarity they show on 
a certain range of the probe length justifies the (str ictly 
speaking improper) use of this term to qualify them fl2l |. 

The description common to all these systems relies on 
the diffusion equation 



d t P{t,r) = DAP{t,r) 



(1) 



written for the probability distribution P (£, r) and 
completed by the conditions imposed on the gradient 
^ P (t, r £ dS) at the boundaries OS of the area S occu- 
pied by the structure at a given moment of time t. In the 
quasistatic regime, when a typical time of boundary vari- 
ations always remains much longer than SD -1 , the prob- 
ability satisfy the Laplace equation AP (t, r) = with 
the boundary conditions [13| evolving as the structure 
grows. This regime allows an elegant description based 
on conformal mapping, which accounts for the structure 
dynamics 0, [l4[ • 

An alternative description relies on random walk mod- 
els, including the widely known Diffusion Limited Ag- 
gregation (DLA) model (l5| , which yields a good visual 
agreement with experimentally observed morphologies 
[8]. The models are based on a probabilistic dynamical 
process 



(2) 



according to which the i-th particle located at a point 
~r^i jumps into the point 7^ within the step t — >• (t + At) 
with the probability p (7^, 7^) = p (7^, 7^). This pro- 
cess is assumed to last until the particle touches another 
particle and stops, thus contributing to an immobile ag- 
gregated dendrite. A link between those two approaches 
was identified in [l6j], giving rise to hybrid methods such 



as the iteration of random conformal mappings 
which generate DLA-type objects. 

Still the results of experiments with atomic clusters 
deposited at a surface and aggregated to dendritic struc- 
tures in the course of difussion, are just in a qualitative 
agreement with the DLA simulations (20|, whereas care- 
ful inspection reveals discrepancies in compactness of the 
experimental and the DLA-induced structure, thus indi- 
cating that some physical ingredients are missing in the 
model. In the present work, we assume that the missing 
element is a Laplacian field induced by the aggregated 
clusters and acting onto the moving ones. The Laplacian 
field can be of different origins, and one of the options is 
a surface stress field, which results from the mismatch of 
the cristalline structures of the cluster and the substrate 
matherial and yields mutual attraction of the deposited 
clusters. In this new setting, the actual motion of clusters 
is no longer a free but a biased random walk with asymet- 
ric jumping probabilities p(7^,7^-,t) ^ p (7^-, 7^, £), 
which on average can be described by the well-known 
Smoluchowski equation for the probability distribution 

d t P (t, r) = (7 (t, r) • ^) P (t, f) +D(t,r) AP (t, r) , 

where the drift force term f (£, r) and the diffusion co- 
efficient D (t, r) are respectively the first and second mo- 
ments of the transition probabilities p. 

The numerical simulation of the resulting "drifted" 
DLA model yields pictures that are in a better agree- 
ment with experimental results. One can go beyond this 
subjective visual estimate and employ the fractal dimen- 
sion as a quantitative characteristic of the dendritic forms 
considered which allows for a meaningful and objective 
comparison between the model and the experiments. We 
show that, even though the Hausdorff- Minkowski dimen- 
sion may be used for our purpose, the so-called external 
fractal dimension is much more adapted due to its great 
sensitivity to the intensity of the stress field. Confronting 
our model to experimental results thanks to this tool, we 
find excellent agreement for antimony-cluster structures. 
Our model, however, suggests that a repulsive force is at 
work in the case of silver clusters, which is incompatible 
with an elastic field. 
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Figure 1: The square lattice in the eight- connectivity model: 
each site is characterized by its position vector r and con- 
nected to its eight neighbours (r + ~et) i=1 8 - 

The paper is structured as follows. In Sec. II, we 
present the drifted DLA model in detail, provide simula- 
tions and discuss the influence of the stress field on the 
morphologies obtained as well as their visual agreement 
with experimentally observed structures. Sec. Ill is de- 
voted to the characterization of the ramified structures 
obtained, either experimentally or numerically, through 
their fractal dimension. In Sec. IV, this quantitative 
tool is used to discuss the agreement between our model 
and some experimental results published in the litera- 
ture. We conclude in Sec. V and suggest some possible 
perspectives of our work. 

II. THE BIASED RANDOM WALK MODEL 
AND THE EMERGING SHAPES OF CLUSTER 
AGGREGATIONS 

We start by presenting some technical details of the 
numerical aproach employed for the cluster aggregation 
modelling and the results of simulations follow. The sur- 
face of the substrate is considered as a two-dimensional 
square lattice of step a. Each site can store one cluster 
at most, and its position is characterized by a vector r 
(see Fig. [Q). 

At the beginning of a deposition experiment, a first 
cluster is dropped onto the surface and diffuses freely on 
the substrate. At each time step of its random walk, 
this cluster can thus jump from the site r it occupies to 
one of its eight closest neighbours of respective positions 

I 



Note that in this model, the clusters are successively 
dropped onto the surface, as it is usually assumed in 
the DLA approach. Though in a real experiment many 



(r + ei), with the same (isotropic) probability p iso — ~. 
At some point of its diffusive exploration of the surface, 
this cluster reaches a default of the substrate where it 
stops. The corresponding site is arbitrarily chosen as 
the origin O of the coordinates - its position vector is 
therefore s\ = 0. This first cluster is now ready to play 
the role of a seed for a ramified fractal-like structure. 

A second cluster is then dropped onto the surface at 
a random site - not "too far" from the origin, otherwise 
it would take it too long to reach the seed by diffusion, 
and diffuses on the substrate now subject to the force 
implied by the interaction with the first cluster fixed at 
the origin. Assuming this force is proportional to the 
elastic deformation of the substrate and therefore scales 
like 1/r according to Hooke's law [21], it can be viewed 
as deriving from the inter-cluster interaction potential 
v (r) = aln(r/ro), where the length parameter ro and 
the energy parameter a characterize the potential range 
and intensity, respectively. The relevant force F (r) = 
— VfW (r) = —af/r 2 acting on the moving cluster affects 
the transition probabilities which now take the form 

P^r + e^l-^F^.e), (4) 

with j3 = l//c#T, as it follows from the Smoluchowski 
equation. For this expression to be valid, the anisotropic 
term —j3F\ (r) • must remain small compared to the 

isotropic probability pi so = |, which ensures that the 
cluster motion never reaches the balistic regime but keeps 
its diffusive character. The second cluster stops its mo- 
tion at a point S2 , as soon as it "touches" the seed located 
at the origin, or in other terms, once it occupies one of 
the eight sites closest to the origin. 

A third cluster is then dropped onto the substrate at 
a random starting position - which, again, should not be 
chosen too far from the two-cluster island already formed. 
This cluster also diffuses on the surface of the substrate. 
Its diffusive motion is however biased by an elastic force 
due to the two-cluster island located around the origin. 
This force is assumed to derive from the sum V (r) = 
cdn (||r — /ro) + aln (||r — S2II A*o) of the potentials 
induced by the first two clusters. The jump probability 
p (r, r + ei) is still given by Eq.(|2]) but now with F (r) = 
— V?V (r). The biased diffusion stops when the third 
cluster visits for the first time one of the neighboring 
sites of the first two aggregated clusters, of position 53. 
The general case of the A^ th deposited cluster follows by 
analogy. 

I 

clusters diffuse over the surface simultaneously, the in- 
teraction among them is negligible as compared to the 
potential induced by the clusters already agregated in 
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Figure 2: Simulations of our ramified structures obtained in our drifted DLA model by the deposition 6400 clusters for different 
values of the dimensionless parameter k: (A)k = 3,3x 10" 3 , (B) k = 1 x 10" 3 , (C)k = 3,3x 10" 4 , (D)k=1x 10" 4 , (E) 
k = 3, 3 x 10" 5 , (F) k = corresponding to the standard DLA model, (G) k = -1 X 10~ 5 and (H) k=-3,3x 10~ 5 . 



a ramified structure. Also note that, in our model, we 
never make use of any specific property of the interaction 
mediated by elastic forces. The nature of the interaction 
might therefore be of a different origin, and the only re- 
quirement is that the interaction potential satisfies 

AV (r) = ac (f) , (5) 

where c (f) = 1 for sites occupied by aggregated clusters, 
and c(r) = otherwise. As a consequence, it is also 
possible to make simulations for repulsive potentials, i.e. 
with a < 0, which do not correspond to elastic deforma- 
tions. 

In Fig. [2] we present a few simulated pictures obtained 
through the aggregation of N = 6400 clusters whose dif- 
fusion was biased by the potential Eq.(|5j) for different 
values of the dimensionless parameter k = a/3. As could 
be expected, for attractive forces n > (Fig. [2]A-E), 
the structures have more compact of the morphologies 
than those obtained in the DLA model k = (Fig. [2]F). 
In the repulsive regime k < (Fig. [2]G,H), the "arms" 
of the ramified structures, on the contrary tend to repel 
each other, thus giving rise to less compact structures. 



III. CHARACTERIZATION METHOD, 
DISCUSSION AND RESULTS 

The fractal dimension appears as the most natural 
mathematical tool to characterize ramified structures 
such as the dendritic islands formed by aggregated atomic 
clusters, and was intensively used in previous works jgj. 
There, however, exist different operational definitions of 
this quantity in the litterature which are not equivalent. 
For instance, the methods based on the radius of gyration 



The transition from sparse to compact objects resembles 
one that has already been observed and modeled in the 
case of electro-deposition experiments |22j], although 
not for the case of self-consistent interaction but for a 
fixed external potential. 

In Fig. [3] for comparison we show dendritic islands 
experimentally obtained through depositing (A) silver 
Agioo [sLU^i an d (B) antimony Sb$oo atomic clusters fl(jl |. 
At first glance, [3] B (antimony clusters) seem to be well 
described by the same kind of structures as in Fig. [2]D, 
corresponding to our drifted DLA model, with a moder- 
ate attractive force, whereas [3] A (silver clusters) seem to 
correspond to none of the structures displayed in Fig. [2j 
i.e. it is satisfactorily described by neither the DLA (Fig. 
[2]F) nore the drifted DLA models (Fig. [2]A-E,G,H). One, 
however, needs a more formal tool for quantification of 
this observation, that would enable one to say how much 
the simulated structures match the experimental data. 
In the next section, we introduce such quantities, the 
regular and the external fractal dimensions. The latter 
presents the advantage of being particularly sensitive to 
the elastic force at the core of our model. 



I 

[iH or the auto-correlation function [l5| as well as mul- 
tifractal analysis techniques (25| appear very sensitive to 
the specific features of the structure under consideration. 
For different DLA simulated morphologies, the fractal di- 
mensions calculated through these methods already show 
a very large dispersion, typically ranging from 1.67 to 
1.72. Because of this dispersion, these methods cannot 
be expected to precisely distinguish simulations obtained 
through different types of models. 

What we are looking for is, on the contrary, a quan- 
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tity able to universally characterize all the morphologies 
simulated by a given model. The fractal dimension Dh 
calculated according to the HausdorfT- Minkowski defini- 
tion precisely has this property. As shown in (26|, it 
takes the value Dh = 1.60 for any DLA simulated struc- 
ture with the negligible dispersion ±0.01, and we nu- 
merically checked the same property for the drifted DLA 
model. To compute Dh one covers the object with disks 
of radius r, and measures the area 7#(r) of the object 
thus formed. The number Njj(r) of these disks is then 
simply given by Nn(r) = • The fractal dimension 
Dh is, by definition, determined by the expected asymp- 
totic behaviour Nn(r) ~ t ~ Dh when r — > 0. Practi- 
cally, the size of the probe disks should, however, re- 
main larger than the typical size r c of a cluster in order 
to avoid spurious and meaningless discretization effects. 
It must also not be too large, that is, typically, smaller 
than the radius R of the whole aggregate - when r ex- 
ceeds R only one disk is necessary to cover the struc- 
ture and therefore A/#(r) = 1. Finally, the quantity 
— Dh = lim ln ^( r ) appears as the slope of the ap- 

proximately linear function ln7V#(r) = /(lnr) on the 
intermediate range [r Cl R] (see FigH|). 

The lower curve in Fig. [5] shows the fractal dimen- 
sion Dh of simulated structures as a function of k. As 
expected, the dimension Dh grows with the intensity of 
the force, which is consistent with the visible increasing 
compactness of the forms numerically generated. The 
variations of Dh, however, remain very moderate: quan- 
tifying the effects of the elastic field through Dh will 
therefore not constitute a very sensitive method. 

To remedy this problem, we consider an aternative pos- 
sible characteristic of a ramified fractal-like morphology, 
the so-called external fractal dimension De, which turns 
out to be much more sensitive to the variations of k. To 
our knowledge, this fractal dimension has never been ap- 
plied in the context of cluster aggregated structures. By 
definition, De is the fractal dimension of the periphery 
of the object. To be more explicit, one computes De by 



measuring the area Ie(?) swept by a disk of radius r, 
sliding along the border of the structure. Intuitively, in a 
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Figure 4: The number of disks Nh (r) covering a ramified sim- 
ulated structure as a function of the probe disks' radius r, in 
logarithmic scales. Self-similarity of the fractal-like dendritic 
structure can be checked on the range [r c , R]. 



ramified object, as r decreases, more and more structures 
appear. This results in a non trivial behaviour of 
and therefore of Ne(t) = ^fpr 5 the number of disks nec- 
essary for covering the whole area i#(r). For a rigorously 
fractal object N E {r) oc t~ De . The same scaling is found 
in morphologies showing self-similarity on a given range 
of radii, as the structures generated by our simulations. 

The behaviours of De and Dh can be compared in 
the inset of Fig. In particular, it is quite apparent 
that De is much more sensitive to k than Dh and is 
therefore a better suited quantity for characterizing the 
modifications imposed to the aggregated structures by 
the elastic field. 
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Figure 5: The inset presents the behaviours of Dh (lower curve) and De (upper curve) of the morphologies obtained by the 
numerical simulation of our drifted DLA model, as functions of the parameter arcsinh (k/ko) characterizing the intensity of the 
elastic field induced by the growing island (k = 0.17 x 1(T 4 ). The main plot is a zoom of the inset on the area contained 
in the dotted box. Two couples of points have been reported, corresponding to the values of Dh and De, extracted from a 
set of experimental pictures obtained by the deposition of silver Agi 00 [9, 23] and antimony S^soo atomic clusters [T(| onto 
graphite HOPG. The value of the parameter k is adjusted so as to minimize \D H 1 (k) — D e ^ p \ + \D t ^ (k) — D e E Kp \ in both cases. 
The points corresponding to antimony can be reported on the theoretical plot with excellent agreement and correspond to 
ftsb — 5 x 10 -6 (small attractive potential), while our model seems to fail to correctly describe the experimental values for 
silver. 



IV. COMPARISON WITH EXPERIMENTS 



Let us now confront our model to some experimen- 
tal data available in the literature. We first consider 
experimental pictures of dendritic structures obtained 
through depositing Sbsoo antimony atomic clusters on 
a graphite HOPG substrate [loj], such as Fig. EflB. Fol- 
lowing the same procedures as described above, we ex- 
tract the Hausdorff and external fractal dimensions of 
the set of ~ 100 pictures we have, and get the aver- 
aged values D e * p (Sb) ~ 1.622 ± 0.022 and D E P (Sb) ~ 
1.876 ± 0.080. If our model correctly describes the ag- 
gregation of antimony clusters, there must exist a spe- 
cific value of the parameter k, denoted by ^sb, such that 
drifted-DLA simulations run with this specific parame- 
ter ftsb yield values for Dh and De which coincide with 
the experimental ones. In other terms, there should ex- 
ist a value ^sb such that D e ^ p (Sb) ~ D H m (ftsb) and 
D G fl P (Sb) ~ D E m (ftsb) are simultaneously checked - 
here, D H m E (k) denotes the HausdorfT/external fractal di- 
mension of the structures obtained by drifted-DLA sim- 
ulations with the parameter k. To identify such a value, 
one simply tries to numerically minimize to zero the 
quantity \D H m («) - D e * p (Sb) | + \D E m (k) - D e * p (Sb)|. 
It turns out that the minimization is successful and yields 



the value ^sb — 5 x 10 -6 . The two experimental values 
D e ^ E (Sb) are reported on Fig. [5j the agreement is excel- 
lent indeed between model and experiment, which sug- 
gests that the physics of aggregation of antimony clusters 
is actually correctly described by our model and involves 
a small attractive elastic potential. It is also interesting 
to note that ^sb — 5 x 10 -6 is precisely the value cor- 
responding to the picture Figj2p: this means that the 
experimental picture Fig|3j3 is indeed better described 
by a structure of the kind of Figj2p than by a pure DLA 
morphology Figj2f\ The visual impression we had at the 
end of Sec. II is therefore completely confirmed by our 
quantitative analysis. 

We now turn to pictures obtained by depositing Agioo 
silver atomic clusters d, (23[. The experimental val- 
ues we extract are D e * p (Ag) ~ 1.528 ± 0.028 and 
D e E p (Ag) ~ 1.677 ±0.138. Note that the error bars, 
especially on D e E p (Ag), are too important to make 
any firm conclusion. This dispersion follows from the 
small size of the dendrites observed on experimental 
pictures, which substantially truncates the range [r Cl B\ 
on which self-similarity can be observed. The fractal 
dimension therefore cannot be extracted with a suffi- 
ciently high fidelity. If, however, we apply the same 
procedure as for antimony, we observe that there ex- 
ists no value of k minimizing to zero the quantity 
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\D°™ («) - D^ xp (Ag) | + |D| m (k) - D e * p (Ag) | . The ex- 
perimental values D e ^ p E (Ag) reported on Fig. [5] corre- 
spond to the value ^A g — — 10 -4 of the parameter k 
providing the best compromise, though clearly far from 
being satisfactory. Note that ^Ag < does not corre- 
spond to an elastic (attractive) field, the physics of ag- 
gregation therefore does not seem governed by the same 
mechanisms as for Sb. More investigations need be taken 
to better understand how to account for this specific be- 
haviour. 



V. CONCLUSION 

We have presented a new class of biased diffusion mod- 
els for the aggregation of atomic clusters deposited on a 
plane substrate. This approach can be viewed as a devel- 
opment generalizing the traditionnally used DLA model, 
that we completed by introducing an elastic force orig- 
inating from the substrate deformation by the growing 
structure. The diffusion equation which governs the clus- 
ter density dynamics is therefore transformed into the 
Smoluchowski equation. Numerical simulations of our 
model yield ramified structures, whose compactness is 
seen to increase with the intensity of the deformation 
force. 

To quantitatively compare the results of our simula- 



tions with experimental data, we investigated Hausdorff 
and external fractal dimensions, and showed that the lat- 
ter is better suited for characterizing the modifications of 
the structures obtained in our model due to the elastic 
field. 

Finally, we compared our simulations to already pub- 
lished experimental results, on silver and antimony clus- 
ters. Whereas our model describes antimony-cluster- 
aggregated structures very well, it is not adapted to 
silver. This suggests that the dynamics of aggregation 
is dominated by different processes, depending on the 
atomic nature of the clusters considered. Identifying 
what these physical processes are, and their exact influ- 
ence on the form of the morphologies of the correspond- 
ing ramified structures is an interesting perspective of 
the present work. Another possible investigation would 
consist to apply the same kind of analysis, involving the 
fractal dimension, to physically completely different do- 
mains, such as electro-deposition where DLA approach 
clearly fails [22]. 
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